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Abstract. The compatibility of the semiclassical quantization of area-preserving 
maps with some exact identities which follow from the unitarity of the quantum 
evolution operator is discussed. The quantum identities involve relations between 
traces of powers of the evolution operator. For classically integrable maps, the 
semiclassical approximation is shown to be compatible with the trace identities. 
This is done by the identification of stationary phase manifolds which give the main 
contributions to the result. The same technique is not applicable for chaotic maps, 
and the compatibility of the semiclassical theory in this case remains unsettled. The 
compatibility of the semiclassical quantization with the trace identities demonstrates 
the crucial importance of non-diagonal contributions. 
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Introduction 



This paper focuses on quantum maps which are represented by unitary evolution 
operators on a Hilbert space of a finite dimension. The quantum map propagates any 
initial state in the Hilbert space by 

V't+i = u^t = c/*+Vo • (1) 

It is assumed that there exists an underlying area preserving classical map which acts 
on a compact phase space M., which can be considered as the counterpart of the 
quantum map within the semiclassical approximation. That is, 

(i) The semiclassical quantization of the classical map (see section ||) provides an 
approximation to the exact quantum map. The semiclassical approximant is unitary 
f/sci(C4ci)^ = I and satisfies the composition rule Ul^i = ^IcPlcX within the accuracy 
margin of the semiclassical approximation. 

(ii) The dimension of the Hilbert space, M, is related to the classical phase space 
volume |A^| by 



M = 



\M\ 



_{2TTh)f_ 

where [•] stands for the integer part, and / is the number of classical freedoms. 



(2) 
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Consider a quantum map U and introduce the notation t„ = tr[/". The fact 
that U is unitary imposes various relations amongst the t„ which should be satisfied 
identically. The purpose of this paper is to study to what extent the semiclassical 
approximation is compatible with a certain class of identities, which, to the best of 
my knowledge, was not examined in this context till now. To give an idea, the most 
simple version of the identities to be considered is 

lim 



ini e ^ t* tn+iyC = 



for arbitrary integers no and v. These identities will be proved and discussed in detail 
in section g . However, before doing this, it is instructive to review a few other 
identities, which involve the tn, and which were used in past investigations of the 
semiclassical approximation ^, ^, ^. 

The first class of identities can be derived by studying the properties of the 
characteristic polynomial 



M 

{z) = det(/ -zU)^Y. ■ 



Since p^{U^) = 0, 

M 

«0 = 1 = - E ^rniU^T ■ (4) 
m— 1 

Multiplying by [/" {n > M) and taking the trace one gets 

M 
m— 1 

By successive applications of the above relation, all the traces t„ with n > M can be 
expressed in terms of the traces of the M lowest powers. 

An important consequence of the unitarity of U is the inversive symmetry of the 
coefficients am, 

flm = e'^alj_^ , (6) 

where det{—U) = e*®. One can utilize the inversive symmetry to obtain identities 
between the tn by invoking Newton's identities. They relate the traces tn and the 
coefficients of the secular polynomial am- 

tm + E aktm-k I ■ (7) 
fe=l / 

Since am = for m > M, the <„ for all n > M depend linearly on the lower m < M 
traces, which is consistent with our previous observation. Successive iterations yield 
explicit expressions for the coefficients am in terms of the t„, and one can substitute 
them in (H) or in (^. 

The significance of such relations in the semiclassical context is due to the fact that 
the tn are expressed semiclassically as sums over n-periodic orbits of the classical map. 
Thus, the compatibility with the exact identities implies that there exist identities 
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relating sums over periodic orbits of different periods which are satisfied within the 
margins of the semiclassical accuracy. The resulting identities between the tn are 
getting complicated as n and M increase, and therefore their compatibility with the 
semiclassical approximation was seldom tested [|j . 

Another set of identities which will be shown to be closely related to the present 
work were introduced by M.V. Berry in He considered the spectral density of a 
quantum hamiltonian 

oo 

d{E) = ^ ^(i? - En) , (8) 
Ji=l 

and its smooth approximant 

d^E) ^ySeiE- En) with 6e{x) = - ■ (9) 

n— 1 

Assuming that the spectrum has no degeneracies (£"„ ^ E^ when n ^ m), one finds, 
27rlim e dl{E) = d{E) . (10) 

Substituting the semiclassical trace formula in both sides of (|o|) one sees that 
the left hand side is quadratic while the right hand side is linear in the periodic orbit 



amplitudes. Integrating (10) over a sufhcicntly large energy domain, the contributions 
of long orbits to the right hand side can be made arbitrarily small, while their 
contribution to the left hand side will not be smoothed out. This observation led 
Berry to conclude that the long periodic orbits must contain information about the 
short periodic orbits if the semiclassical approximation is compatible with ( p^ ) . This 
information is stored as correlations between actions of periodic orbits, because pairs 
of distinct periodic orbits combine together to give an amplitude of order in 
d'^{E), which, upon multiplying by e reproduce the periodic orbit contributions to 
the oscillatory parts of d{E). This is a highly "non-diagonal" effect, which needs very 
special correlations between the actions to get the correct result. This observation 
shows that the use of identities of this type comes naturally in the context of the 
study of classical action correlations and their effect on the statistics of the quantum 
spectra §, [l0[ 0. 



J. Keating [[7| (see also generahzed ( |lO|) and used it in his studies of the 
spectrum of the Riemann zeros. The "non-diagonal" correlations which are necessary 
to prove the identities are introduced by using the Hardy-Littlewood conjecture on 
the correlations between primes. 

E. Bogomolny tested ([l0| ) for the spectrum of a rectangular billiard with 
periodic boundary conditions. By considering carefully the stationary phase manifolds 
in the sums over periodic orbits which arise from the left hand side, he was able to 
perform the summations and to demonstrate the compatibility of the semiclassical 
approximation with (^^. The same methods were later used to analyze integrable 
systems in general. The present work was inspired by a seminar given by Professor 
Bogomolny, and it can be considered as an extension of Bogomolny's ideas to quantum 
maps and their semiclassical approximation. 

This manuscript is arranged in the following way. The trace identities which will 
be the main tool in the present analysis will be derived in the next section. The 
semiclassical quantization of area preserving maps will be reviewed in section ^ and 
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the compatibility of the semiclassical approximation wiU be demonstrated in section 
^ for integrable maps. The close relation between the compatibility problem and the 
correlations between actions of periodic orbits will be addressed at the end of this 
section. The difficulties encountered in attempting to use the same methods to study 
the compatibility in the case of hyperbolic maps are discussed as well. By applying 
the trace identities to certain quantum graphs, it is possible to derive combinatorial 
identities which involve Krawtchouk polynomials These polynomial are 

important building blocks in the theory of error correcting codes p5| ] , and the identities 
might be of use in this branch of mathematics. Since this application is not directly 
connected to the issue of the semiclassical compatibility, it will be presented in the 
appendix. 

2. Trace identities 

Consider a unitary matrix U of dimension M. Its spectrum consists of M points on 
the unit circle {e*^™, m = 1, • • • , M}, where the eigenphases are real and are assumed 
to be distinct. Recalling the notation t„ = trC/", the following identities hold for 
arbitrary integers Uq and v : 

oo 

lim e J2 tntn+^e-"" = U ■ (11) 

n— no 

To prove these relations one substitutes t„ = X]m=i e*"^™ and after summing the 
geometric series, one uses 

e (I for i0,n - e,n') = \ _ . 

e'™o 1 - c'(»'"-e™')-e - \ for (0,„ - e,n') ^ / " ■ 

The condition that there are no degeneracies in the spectrum of U is used to justify 
the rightmost equality in (|l2|) . 

A few points are worth noticing: 

- The e — *■ limit of the sum weighted by e""*^ can be interpreted as a time average 

lim e E(-)ne-^Jm^^ (13) 

n— no n— no 

- For 1^ — 0, and using ~ M, one gets 

oo ^ 
n=0 

Thus, the time average of ;p'|tnp approaches 1, a result familiar from the study of the 
spectral form- factor for unitary matrices [ p6| . 

- The spectral density of U on the unit circle can be written as 

M oo 

d{e) - V 5(0 - Oi) = lim d,{e) ■ d,{e) = — V t^e-^^^-^l™! . (15) 

l—l m— — oo 



Using (|llD, one can easily derive 



27rlini e di{e) = d{6) , (16) 
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which is the analogue of (|T^) for the spectrum of unitary operators. 

The following identities which involve products of more than two traces can be 
proven with the help of (|l^ ) 

= lim • • • lim ei ■ ■ ■ Ck (17) 

oo C30 k 

■■■ ^ d'"Cfc*(riiH hnk)+i' exp(— ^njEi), 

where the are arbitrary integers. The identities (17) can be also written as 

lim • • • lim ei • • • et (18) 

C30 oo k 

'■' ^riiH — i-nfc^"i+i'i ' ' ' ^"fc+i^fc cxp(— n^ei). 

The equivalence of (17) and (18) can be shown by taking the complex conjugate of 
(18), denoting vi + ■ ■ ■ + Vk = —v and shifting the summation variables ni by Vi. The 
shifts do not affect the result since they can be absorbed in the arbitrary {n^g}. 

The compatibility of the semiclassical approximation for tn with the identities 
(pi], 17) will be investigated in section ^, after the semiclassical approximation for i„ 
will be reviewed in the following section. 

3. Semiclassical quantization of maps 

The quantum unitary operator \J to be investigated, is assumed to be the analogue of 
an area preserving map T acting on a finite phase space domain M. with area \M.\. 
(For the sake of simplicity the maps to be considered act on manifolds with / = 1, 
and will be assumed to have the twist property. The semiclassical treatment can be 
extended to the general case.) The phase space coordinates are denoted by 7 = (g,_p) 
and 7 is mapped to 7 = T(^^\ Area-preserving maps can be derived from a generating 
function (action) $(g,g) 

The explicit mapping function 7 = ^^(7) is obtained by solving the implicit relations 
(|l9|). The twist condition ensures that the implicit equations ( |l9[ ) have a unique 
solution. 

In the case the map is integrable, let / be the invariant momentum (action 
variable) under the action of the map and the canonically conjugate angle variable. 
The domain of the map is the annulus / € [Imim Imax], 4* G [0, 27r), and \M\ = 
2'K^Imax — Imin) ■ In this representation, the generating function must take the form 
i>((/), ^) = $(^ — 0). The explicit map is 

/ = / ; A0 = 0-0 = /(/) (20) 

Here, /(/) (the angular velocity) is the inverse of the generating relation / — $'(A0), 
which gives A(j> in terms of /. The twist condition is fulfilled if $"(A(^) 7^ 0. 

A classical trajectory is obtained by applying the map to an arbitrary initial point 
in M , and the corresponding action is accumulated along the trajectory. 
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The semiclassical quantization of T in the q representation is Wi 



(9|t^l9)sci = 



/ 1 y 


[92$(<7,g-)1 


ei*(9'9) 


yiTihi ) 


dqdq 





(21) 



It can be shown to preserve the composition property and to be unitary within the 
semiclassical approximation. 

The semiclassical approximation for t„ involves the periodic manifolds of the 
classical map. For hyperbolic maps ||, 

pt^„ |det(/-r;)|2 

The semiclassical approximation for tn involves the set of n-periodic orbits Vn which 
are repetitions of primitive periodic orbits of T ^ with periods rip which are divisors of 
n, so that n = n^r. The monodromy matrix is denoted by Tp. Each periodic orbit 
contribution is endowed with a phase which is the action summed along the periodic 
orbit, 

Tip 

*p = (with g„^+i = gi), (23) 

and of the Maslov contribution —Vyix jl. 

For integrable maps, the action-angle variables (/, will be used, where J is the 
classical invariant. In the quantum picture, / is quantized to integer multiples of h, so 
that Ij — jh and 1 < j < M. The matrix U is diagonal in the action representation. 
The semiclassical approximation for the eigenphases can be carried out directly, 

where /(/) is the angular frequency (|20|). Thus, 

M 

[i„]sci = 5]e™[^^l- , (25) 
i=i 

This semiclassical expression is not of the desired form, because it does not express 
t„ in terms of the periodic orbits. However, performing the j sum using Poisson 
summation, one gets, 



m— 1 



Now, tn is expressed as a sum over the periodic manifolds of period n and winding 
number m. They occur at values of I for which the angular frequency is rational 

f{In,m) = 277^. 

The expressions for i„ in terms of periodic manifolds in the cases of classically 
integrable and classically chaotic maps are the necessary building blocks for the 
discussions which follows. 
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4. Compatibility of the semiclassical approximation 

The compatibility of the semiclassical approximation for the t„ with the trace identities 
will be now shown for integrable maps. Turning first to (O), the semiclassical 
expression (Qm for t„ is substituted in the right hand side of (nW resulting in 



n— no 



2 

n 



(27) 



(n{n + ly))'^ 



m=l m' = l 



$"(2^_)$"(2^— 

n n - 



,^[("+!^)*(2'r^)-r!*(27r^)] 



The sum above runs over the periodic manifolds (tori) of the map. It is important 
to respect the integer character of n,m,m', since only for integer values, the classical 
orbits are periodic. If we turn these sums to integrals by e.g. Poisson summation, we 
would loose this feature. To get a finite contribution for (27) we must collect all the 
terms which contribute coherently to the sum. Summing over these terms (weighted 
by e~'^") should provide a pole at e = 0, so that the final multiplication by e will yield 
the residue at the pole. Inspecting (27) we immediately notice that, for example, the 
terms for which = = a{i^) (where a(j^) depends only on v) yield a contribution 
of the desired nature, because the net phase i^$(27ra(i/)) is common to all the summed 
terms. To find these contributions in a consistent way, we shall identify them as the 
points where the first variation of the phase of the summand vanishes. 



5 {phase} — Sr 



$(27r- 



_)_$(27r-) 
ly n 



5m' 



27r 



27r$'(27r. 



m! 



-$'(27r 



— $'(27r— 

n n 



Sm 



27r$'(27r 



TO, 



The phase is stationary with respect to variations in n when 



TO 

n 



(28) 



(29) 



One should consider only the solutions in the range m < n and m' < n+v, consistently 
with the range of the sums over to and to'. The solution of (29) in integers is 



= Nv 
= Nk 



with 
with 
^ l)fc . 



N = No, 
k= I,-- 



■ ■ ,oo 



(30) 



The arbitrary integer hq is replaced by another arbitrary constant Nq which fixes the 
lower limit of the n sum. This solution is the general solution for the generic cases. 
One can always invent maps for which other stationary points exist. For each value 
of k, the points of stationary phase form a grid. Near each grid point, the summation 
variables will be replaced by local variables 



n = Niy + 6n : m — Nk + 6m 



{N + l)k + Sm' , 



(31) 



8 



with 

l^nl < \ ■ < \ ■ \6m'\ < ^ . (32) 

The range of variation of the local variables is chosen such that each point in the 
original sum will be counted once. The contribution of the domain about each grid 
point will be computed by using the stationary phase approximation. The phase 
cannot be made stationary with respect to independent variations of m and m' because 
there is no guarantee that $'(27r-2^) ~ and <i>'(27r— ) = have solutions when n , 
m and m' are integers. However, when (^9|) is satisfied, the phase is stationary on the 
manifold Sm = 5m' . Thus, the sum over n, m, m! in (27) is replaced by 

OO n n+u p oo \ 2 2 2 I 

EEE-EE E E E ^^-m. (33) 

n=no m=l m' = l k=l N=No [5n=-f 5m=-^Sm'=-^ J 

Where the curly brackets on the right enclose the contribution of the vicinity of a 
single grid point, restricted to the line 6m = Sm' . On this line the summand is 
constant and therefore the summation amounts to multiplication by k. The Sn sum 
can be computed by considering the quadratic approximation to the phase near the 
stationary points and retaining the leading term in the result. It is determined by the 
curvature of the phase at the point where it is stationary 



9^{phase} 



The amplitude of the result depends on N, but when it is substituted in (27) it exactly 
cancels the N dependent coefficient. The phase of the Sn sum is i/$(27r^) which is also 
independent of N. Thus, the resulting terms of the N sum depend on N only through 
g-ei^AT^ Summing and taking the limit e ^ results in a factor Collecting all the 
factors, one finds that (27) is approximated by 



lim , J2 K.'"+-l.cl<''"' 



(I) ^ 

^ ^ k=l 



$"(2^-) 

V 



e^'^*(2-#) = [i,]^^, . (35) 



This completes the proof that the trace identities ( |ll| ) are compatible with the 
semiclassical approximation. 

The more complex relations (17) or the equivalent (18) can be checked using the 
same technique. As an example, the compatibility of the identity 

lim £162 / ^ d+n2*"i+''i*"2+i'2e "^'^ "^'^^ = ^1/1+1^2 (36) 

ei,e2— U ^ — ' 
"1 ,"2 

will be demonstrated in some detail. 

Substituting in (^6|) the semiclassical approximation ( p^ ) for one has to consider 
the sum 



fi ) ^ \ni + n2 Til + I'l n2 + 1^2 , 



rii+i/i 712 + 1^2 ni+n2 

<^"{2^ )$"(27r )$"(27r- 



E E E 



X 

^ ^ 1 I ' ni + i^i ' ' n2 + V2' ' ni + n2 ' 

mi — 1 m2 — 1 mi2 — 1 
X ei[("i+'^i)*(2'^;rfTT^) + ("2+!^2)*(2x^;|i2;^)-(rii+n2)*(27r^j^^)] ^-^^^ 



The phase of (37) can be made stationary with respect to variations of ni and n2 
under the following conditions: 

- $12 + $1 + ^$'i2 = 38 

ni + 712 ni+vi 

Where the short-hand notation $i = $(27r;j^2ii_) and $i = 27r$'(27r;j^^), etc was 
used. The conditions (38) are satisfied by solutions in integers of: 

= < 1 and = < 1 . (39) 

ni +712 «i + t^i ni +n2 n2 + V2 

The general solutions depend on three integers A^i, N2 and k so that 
ni = {Ni-l)vi+ N1V2 with Ni^Na,---,oo 

7l2 = A^2'^l + (^2 - 1)'^2 with A^2 = ^0, ' ' ' , OO 

TO12 = [Ni +N2- l)k with 1 < fc < I'l + 1^2 (40) 

777,1 = -^1 k 
7712 = -^2 k 

Again, for each value of k the points of stationary phase form a grid and one computes 
seperately the contribution from the vicinity of each grid point. The summation 
volume about each grid point is of size 

\6ni,Sn2\ < ^^^^^ ; |(5mi, 57772 1 < ; |(5777 i2| < fc (41) 

The phase cannot be made stationary with respect to independent variations of 
Till, 7772 and 77112. However, like in the previous case, the phase is constant for 

5mi2 = 6mi + 6'm2 (42) 

The sums in (37) are rewritten in terms of the local variables for each grid point, 
and the curly brackets encloses the sums over individual domains (^). 

00 00 ni + i/i n2 + i^2 ^^l+^^2 

E E E E E - (43) 

m— no n2— no mi — 1 ni2 — l mi2 — 1 

E E E ^ E EE E'^'5™i2.*'"i+'5™4 

fc=l Ni=No N2 = No {Sni Sn2 Smi Sm2 ) 

The last Kronecker S is due to the restriction (p2|). Since the summands in the 
Smi,Sm2 sums are constant, they contribute a factor of k^. The S'ni,Sn2 sum is 
performed again by expanding the phase to second order, and retaining the leading 
term in the result. The determinant of second derivatives at the point of stationary 
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phase is 



det 



( 



52 



{phase} 



) 

4 



(44) 




drLidrij 
k 



) 



Vi + V2 



(rii + n2){ii\ + vi){n2 + V2) 



where n\ and 712 take the values (40). Collecting all the factors, and performing the 
A'^i and N2 sums while taking the limit ei, £2 — ^ 0, one remains with the sum over fc 
which can be easily identified as [^i/i+iyalsd- 

This method can be extended to identities involving higher powers. The procedure 
becomes much more cumbersome and will not be reproduced here. 

The essence of the derivations outlined above is that the phase of the summands is 
constant on an infinite grid of integers. Only when all of them are summed together, 
they provide the necessary singularity which is cancelled against the e factors and give 
the correct answer. 

The compatibility of the trace identities with the semiclassical approximation 
demonstrate the importance of "non-diagonal" correlations which are essentially due 
to the repetitive nature of the distribution of periodic orbits for intcgrable maps. 
As a matter of fact, had one applied the standard diagonal approximation where 
repetitions are neglected, to the sums (27), (37) one would get a vanishing result when 
V or {y\ + 1^2) ^ ■ The condition ( p9|) picks up (non-diagonal) pairs of n-periodic and 
n + v periodic manifolds which coincide, since the action variable In,m is the same. 
This feature is effective in intcgrable systems since the periodic manifolds are specified 
completely by integers (the period n and the winding number m),and it is responsible 
for the compatibility of the trace identity with the semiclassical approximation. The 
condition (40) expresses a similar coincidence of three periodic manifolds. In chaotic 
systems, repetitions do exist but they play a much less important role, because apart 
from the period n there exists no other integer which would replace m to specify the 
periodic orbits. This is why the methods described in the present section fails for the 
chaotic case. The substitution of the expressions (j2^) for systems which are chaotic in 
the classical limit, and attempts to identify the classical correlations which are implied 
by the trace identities failed, so far, to give a definite answer. 

5. Discussion 

In the previous section, the compatibility of the trace identities with the semiclassical 
quantization of intcgrable maps was studied from the point of view of periodic 
orbit theory. It can be also proven very simply by observing that the semiclassical 
expressions (^5|) for the t„ has the formal structure as a trace of a unitary matrix, 
and therefore the trace identities which are based on nothing else but on this form, 
are automatically ensured. This result does not detract from the work presented in 
the previous section, because it gave new insight into the interplay between the trace 
identities and the correlations in the spectrum of classical actions. 

In general, the semiclassical approximation for a quantum evolution operator is 
the leading approximant in an asymptotic expansion, in which the time is kept fixed 
and h approaches zero. Taking the opposite order in which h is fixed and the time 
is increased, is not allowed within the semiclassical approximation. One has to bear 
this basic rule in mind when one tests the semiclassical approximation. In the present 
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context, the magnitude of sets the scale of the relevant evolution time (see (]l3|)). 
Hence the proper order of limits would seem to be such that e be fixed while ?i — > 0. 
However, the semiclassical limit is synonymous to M — > oo (^. In this limit the 
spectrum of U becomes dense on the unit circle, and the trace identities do not hold. 
These conflicting demands can be satisfied when e < 0{h-^). This remark does not 
apply to the test we performed for integrable maps, because the formal structure of 
[in]sci as given by is compatible with the trace identities for any values of h and 
n. This is not the case for quantum maps which correspond to classically chaotic 
systems. 

So far, any attempt to assess the compatibility of the trace identities with the 
semiclassical approximation for classically chaotic systems ended in failure. Even for 



maps on graphs, where the periodic orbit expansion of i„ is exact |19| , |20| , 21 , I was 
unable to derive the trace identities within periodic orbit theory. Therefore, the nature 
of the underlying classical correlations which ensure the exact trace identities when 
they are expressed in terms of periodic orbits remains an enigma, and calls for further 
research. In Appendix A, the trace identities are used to derive some combinatorial 
identities, which illustrates an application outside of the domain of periodic orbits 
theory. 
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Appendix A. Identities for Krawtchouk polynomials 

We consider here the quantum map for a 2-star graph which was defined and studied in 
detail in |^ . It consists of a "central" vertex out of which emerge 2 bonds, terminating 
at vertices (with indices j = 1,2) with valencies Vj — 1. The ratio between the bond 
lengths Lj is assumed to be irrational. This simple model is not completely trivial if 
the central vertex scattering matrix is chosen as 

^(o) = [ 75 ^ ) . (Al) 



Neumann boundary conditions are imposed at the two other vertices. The Hilbert 
space is of dimension 2 and the evolution operator is 
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where the diagonal matrix on the left takes care of the free propagation on the bonds 
and the reflections from the vertices j = 1,2. 

One can write an exact expression for t„(fc) = trJ7"(fc) in terms of periodic orbits 
on the 2-star graph which is analogous to (p2|). There are 2"/ri n-periodic orbits. 
However, their lengths can take only n+1 distinct values : L{n, q) — 2(qLi + (n—q)L2), 
with < q < n. Thus, one can write 

n 

i„(fc) =^e^^("'«)'=A(n,g) , (A3) 

9=0 

where A(n, q) is the coherent sum of the amplitudes contributed by the isometric 
orbits with length L{n, q). It can be computed explicitly [ pT| in terms of Krawtchouk 
Polynomials 

I ( 1 for g = or n 

= (-l)"+.yf(pp„_i^_,(g) for 0<g<n. ^^^^ 

and the Krawtchouk polynomials are defined as in |l^ by 

P.,.(.) ^ (^) jZi-ir- Q (^:;) for < < . (A5) 
Substituting ([A^ ) in the trace identity ([III), we get for any integers u and uq 

oo n n-fi/ 

lim e J2 e-"'^^^e2*'=[(P-«)-^i+(''-(P-«))^=^U(n + zy,p)A(n,g) 

n—no q—0 p—Q 

= ^ K)e2''=[''^i+(''-'")^^l . (A6) 

This is valid for any value of the wave number fc. When the lengths L\ and L2 are 
incommensurate, the equality can hold only if the coefficients of the phase factors 
^ikL(v,K) both sides are equal. This implies 

00 n 

lim e ^ e"'" ^ A{n + v,q + K)A{n, q) = A{v, k) , (A7) 

n=no 5=0 

with < K < ly. Restricting to k to the interval 1 < < ~ 1 the identity can be 
expressed in terms of Krawtchouk polynomials exclusively. 



- ( P.-r,_(.) = lim e ^ ^ ^ (^8) 



X 



q=0 



q{q + k) \qj \q + k 



I could not find identities of this kind in the standard books on Krawtchouk 
polynomials. 
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